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Estimating the abundance of organisms is an 
important aspect of ecology. In fact, if we adhere to 
Krebs' (1972) definition of ecology as "the scientific 
study of the interactions that determine the 
distribution and abundance of organisms" it is 
fundamental to the field. Estimates of population size 
form the basis for ecological and conservation studies. 
A multitude of methodologies exist to estimate 
population abundance. These methods differ in their 
suitability for species, the assumptions involved, the 
accuracy of the estimates, and the effort and cost 
needed to perform. The most appropriate technique 
will likely depend on the objectives of the study and a 
balance between the precision and parameters needed 
and the cost and effort of each method. Because of 
differences in methodology, estimates of abundance 
may not be directly comparable among studies. Here, 
we examine the relationships among several common 
population estimation techniques used for butterflies. 

Butterflies are popular study organisms for a variety 
of ecological and evolutionary questions and many 
species are often used as indicators in conservation 
studies (Blair 1999, Brown and Freitas 2000). It is our 
hope that this study will allow more meaningful 
comparisons of population data collected using 
different methods and provide guidance in selecting 
among common techniques. 


MATERIALS AND METHODS 


Study speeies and site. The butterfly Parnassius 
smintheus Doubleday (Papilionidae) is abundant 
within subalpine meadows in the Rocky Mountains, 
although congeners are threatened elsewhere (Kuras 
et al. -2000). The butterflies! host plant, Sedum 
lanceolatum "Torr. (Crassulaceae), occurs in gravelly 
sites above tree-line (Fownes and Roland 2002). 
Parnassius smintheus is univoltine with a flight period 
from mid July to September in our study area. Adults 


nectar on yellow flowered species such as S. 


lanceolatum, Solidago multiradiata (Asteracae), and 
Senecio lugens m e (Matter and Roland 2002). 

Transect surveys and mark-recapture of P 
smintheus were conducted in nine meadows within a 
network of 21 meadows located along Jumpingpound 
Ridge, Alberta, Canada (51? 57'N, 114° 54' W). Each 
meadow was considered as containing a 'population.' 
Meadows are comprised of grasses, sedges, and 
wildflowers, and are bordered on their lower gs by 
forest consisting of Pinus contorta, Abies lasiocarpa, 
and Picea engelmannii. 

Population estimation methods. For transect 
surveys, each observer walked a path through the 
middle (along the longest axis) and around the 
circumference of a meadow tallying the number of P. 
smintheus observed at any distance in front of them. 
As P. smintheus fly more frequently when it is sunny 
(Ross et al., in press), observations were conducted 
during full sun. As a rubric, we stopped walking and 
counting if we could no longer see our shadow. For 
each survey there were between two and four 
observers. Transect surveys were conducted prior to 
mark-recapture, on the same date, to provide 
comparisons. 

For mark-recapture estimates, we captured 
butterflies using hand nets and each newly captured 
butterfly was given a unique 3-letter code on the 
upper surface of each hind wing using a felt-tipped 
pen. For all captures, we recorded the date, time, 
location, and identity mark (Roland et al. 2000, Matter 
and Roland 2002). To equilibrate effort among 
populations, recapture continued until ~75% of 
recaptured buttcrilies had been previously captured 
that day. Populations were sampled from 3-7 times 
between July 27 and August 20, 2001. 

Transect surveys. Transect surveys are perhaps the 
simplest population estimation technique. This 
method assumes, if multiple observers or observations 
are involved, a consistent path or amount of time is 
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used, and that observers have similar ability in 
identification (Pollard 1977, Thomas 1983). To arrive at 
a population estimate for transect surveys, we calculated 
the mean and variance of the number of butterflies 
reported by the observers. 

Number of individuals captured. This was the 
simplest mark-recapture technique. For this estimate 
we tallied the total number of different individuals 
caught during a sampling session. This and all other 
methods Quse ing marked individuals (below) assume 
that marks are not lost, and that marking and handling 
do not affect behavior, survival, or the probability of 
capture. 

Craig's Method. Craig's method is a slightly more 
complicated mark-recapture technique based on the 
frequency of capture during a single sampling session 
(Craig 1953, see also Southwood 1994). Thus, it uses 
some data often discarded by other methods, and may 
not be applicable for small mammals or other organisms 
which are usually caught only once per session or with 
unequal capture probability (Edwards and Eberhardt 
1967, Nixon et al. 1967). Population estimation assumes 
that the frequency of butterflies captured once, twice, 
thrice, etc. follows a Poisson distribution. The number 
of butterflies not caught, the zero term of the 
distribution, is estimated and 'added' to the number of 
individuals caught to arrive at an estimate of population 
size. This method incorporates all the assumptions of 
marking and further assumes that all individuals in the 
population are at equal risk of capture at all times, i.e. 
there is instantaneous re-mixing upon release and no 
handling or marking effects that would effect capture. 
Craig's ‘method also assumes that the population is 
closed, that is there is no birth, death, or migration 
during sampling. Population size was estimated using 
the equation: 


InN -In(N—r)=s/N 


where Ñ is estimated population size, r is the number 
of individuals captured, and s is the total number of 
captures (Craig's method 1, Craig 1953). We solved the 
equation above using the fsolve routine of Maple V. 
Variance of the estimate was calculated as: 


2 
gU TERE 


where A = s/8 (Southwood 1994). Capture probability 
can be estimated as p = r/Ñ. Given S and A, the expected 
number of individuals caught x times can be calculated 
from the Poisson equation: 
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Goodness of fit was evaluated by comparing these 
expected values to the observed values, where 
(observed -  expected)/expected follows a X 
distribution. Evaluation can be made for each class of 
number of captures with one degree of freedom, or 
overall, by summing capture classes with degrees of 
freedom equal to the number of summands. 

Geometric Modcl. Similar to Craig's method, the 
geometric model is also based on the frequency of 
'apture and assumes a closed population. However the 
assumption of equal capture probability is modified and 
the model treats the number of times that an individual 
is captured as a geometric distribution. Population size 
was estimated using the equation: 


m r(s—1) ; 

By e 

and variance of the estimate as: 
Ng/dq' 

where g = (s-r)/(s-1) and 8 = 1- d (Pollard 1977). Note 
that & is used rather than the traditional p to avoid 
confusion with capture probability. Capture probability 
can be estimated as D = r/&. Given Ñ, ĝ and 8, the 
expected number of butterflies caught x times can be 
calculated as: 


E(f,)- N&p^. 


Gooduess of fit can be evaluated. using the same 
methods as for Craig's method (previous section). 

Lincoln-Petersen. As opposed to the previous 
methods, the Lincoln-Petersen method requires 
captures on multiple Occasions, in our case consecutive 
days. This method is based on the assumption that the 
ratio of marked individuals to the total population size 
will equal the proportion of marked individuals in a 
second sample. It assumes the assumptions for marked 
individuals, that populations are closed during and 
between sampling periods, and a constant capture 
probability. The equation: 


was used to estimate population size, where m is the 
number of individuals marked on the first occasion, r is 
the number of recaptures, and n is the total number of 
individuals captured on the second occasion. For 
samples under 20, we used a small sample 
approximation (Baily 1952): 
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Variance of the estimate was calculated as: 
Y: 
pnt 
CoS eee 


r? 


and as: 

> m(nt+l(n—r) 

" (r +1)’ (r+ 2) 
for the small sample approximation (Southwood 1994). 
Estimated capture probability during the recapture 
period for the can be calculated as fj = r/m (Skalski and 
Robson 1992). 

Jolly-Seber. The Jolly-Seber method is similar to 
Lincoln-Petersen, but requires capture on three or 
more occasions. Importantly, this method relaxes the 
assumption of a closed population. Animals may enter 
the population via immigration or birth and leave the 
population via emigration or death. Without additional 
information, estimates can only be made for the 
combined effects of each, that is, total gain and loss to 
the population. The model is stochastic assuming that 
there is a probability that organisms will survive (not die 
or emigrate) from each census period to the next and 
that capture probability may also vary. Survival (6), 
capture probability (5), and population size (Ñ) were 
estimated using the program Jolly. We assumed fully 
parameterized models (time varying capture and 
survival probabilities) unless simpler models with 
constant survival, constant capture probability, or both 
constant did not significantly differ from the full model. 

Analysis. We regression to build 
predictive relationships between the population 
estimation methods. We constructed a separate model 


used linear 


for each pair of methods. Because of non-linearity 


between some estimates, data were log, transformed 
prior to analysis. Standard. diagnostic techniques for 
regression were used including inspection of residuals 
md outliers. Not all — estimation techniques 
could be used for each sample date, e.g. sampling 
needed to be conducted on consecutive days for 
Lincoln-Petersen estimates. Thus, sample size varies 
among the techniques. We considered each population 
estimate to be an independent observation. It should be 
noted that some relationships involve cases where the 
dependent and independent variables are calculated 
using the same data (e.g. Lincoln-Petersen and Jolly- 


Seber both incorporate captures in the estimate of 


population size). In such cases correlations will be 
greater than expected by chance, affecting statistical 
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inference; however, the regression equations describing 
the relationships are still v alid. 

Estimates for a population of known size. To 
estimate the accuracy of the techniques, we released a 
known number of male butterflies (24) into a meadow at 
lower altitude where they had never been observed and 
their host plant docs not occur, but many of their nectar 
flowers do occur. Butterflies were released onto a nectar 
source at varying positions throughout the meadow. 
Sampling began 30 min after release. Butterflies were 
marked and recaptured as in the population surveys. 
Three observers who did not know the number of 
butterflies released, conducted the transect surveys and 
mark-recapture. We conducted one transect survey and 
two mark-recapture sessions separated by one hour for 
this population. We computed population estimates as 
for the natural populations. As there were only two 
capture sessions Jolly-Seber estimates could not be 
calculated. 

RESULTS 

There were significant, positive correlations among 
all the population estimation techniques (Table ii 
Fig.1). Transect surveys produced the lowest estimates, 
while the geometric distribution provided the highest 
estimates of population size. Models for which a test 
could be preformed showed no lack of fit except for 
Craig's method for meadow Z 

For the population of known size (24 butterflies), the 
mean of the three observers! transect counts was 7.3 + 
5.3 (Var.). There were 16 and 14 captures for the first 
and second census, respectively. Craig's estimate for the 
first census was 24.4 + 28.9 and 26.3 + 69.5 for the 
second. The estimates from the geometric distribution 
were 40.0 + 150.0 and 44.3 + 304.2. The Lincoln- 
Petersen index estimated population size as 18.5 + 3.5 
butterflies. 


DISCUSSION 


The significant, positive relationships among the 
population estimation techniques were reassuring. Our 
limited investigation of the accuracy of the tec niques 
shows that transect counts and the number of captures 
underestimate the actual population size. Craig's 
estimates were accurate while the Lincoln-Petersen 
estimate was lower than the actual population size, but 
provided a reasonable estimate. The geometric model 
overestimated population size. This experiment also 
allowed us to test our mode! and illustrate its utility and 
limitations. Note that the prediction of a single value of 

Y and its error for a given X in regression (prediction 
Em differs from, and is greater than the 
distribution of Y (confidence interval) at a particular N 
(Zar 1999, p. 341). As an example, our trausect count of 
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7.3 results in predictions of 26.0 + 4.6 (95% P.I.) for 
Craig's estimate, 12.9 + 3.5 for the number of captures, 
and 32.5 = 3.9 for the Lincoln-Petersen estimate. The 
actual estimates fall within the prediction intervals for 
Craig's estimate and the number of captures, but not for 
the Lincoln-Petersen method. This demonstration 
illustrates both the utility of our model and_ its 
difficulties. For small population sizes it may be difficult 
to obtain a precise estimate. This problem can especially 
be seen by the fact that the intercepts of some 
relationships were significantly different from zero 
(Table 2). For example, a transect count of zero may 
indeed indicate the presence of no butterflies, but could 
result in an estimate of 6.6 based on the Lincoln- 
Petersen estimate. 

Although the equations presented here apply only to 
P. smintheus at our study site, the results do illustrate 
the relative strengths and weaknesses of the various 
techniques. Given the varying reasons for estimating 
population abundance, a variety of methods have been 
and will continue to be used. For butterflies, transect 
surveys are perhaps the easiest and least disruptive 
method of population estimation requiring only the 
ability to identify species on the wing. For some groups 
or assemblages this may be quite difficult, necessitating 
either netting or grouping of species that cannot be 
distinguished. For conspicuous, easily identifiable, 
species transect surveys are an efficient means to 
generate relative estimates provided observability 
(capture probability in mark-recapture terminology) 
does not vary. However, transect surveys do not provide 
accurate estimates of population size, nor do they allow 
for estimation of obscrvability which limits their utility 
for comparison. Transect surveys, as conducted here 
under highly favorable conditions, result in large 
underestimates of population size. This underestimation 
is especially important in determining presence or 
absence. A transect survey producing no butterflies does 
not mean that the species is absent. Accurate 
determination of local absence or extinction always will 
require additional, intensive sampling. 

All other methods investigated here require both the 
capture and marking of individuals which may alter 
behavior (Mallet 1957). In general, capture temporarily 
reduces the propensity of butterflies to fly. Reduced 
flight in turn lowers the capture probability of marked 
individuals relative to unmarked individuals for the 
length of time that the butterflies are affected (Gall 
1955). Frequency of capture methods (Craig's and 
geometric) will be more influenced by a temporary 
change in behavior than other methods. If marked 
individuals temporarily have a lower capture probability 
than unmarked individuals, estimates of population size 
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will be higher than the actual population size (Gall 
1985). For other mark-recapture methods, anv 
temporary handling effect usually will have abated by 
the next census al For P. EL T the effects of 
marking on flight are minimal; capture 
probability is lower for females than for inales violating 
the assumption of equal capture probability for Craig's 
and the Lincoln-Petersen methods (Roland et al. 2000, 
Matter and Roland 2002). We note that our estimates 
for a population of known size used only males, and thus 
should not violate this assumption; however, this bias 
will affect estimates for the populations along 
Jumpingpound Ridge. 

The number of inn iduals captured underestimates 
population size as the capture rate rarely nears 100 
percent. However, assuming marks are not lost, the 
number of individuals captured does provide an 
estimate of the minimum possible population size. For 
the effort of tallying the number of times each 
individual butterfly is captured, Craig's method provides 
fairly accurate estimates of population size at a specific 
time, while thc geometric distribution overestimated 
population size. Iein both frequency of capture 
methods showed good fits to the data despite assuming 
different distributions for capture frequency. In general, 
fits werc better for Craig's method than for the 
geometric distribution. This result contrasts with 
Pollard (1977) who found better fits for the geometric 
distribution than for the Poisson distribution of Craig's 
method in his investigation of three butterfly species. 
Our result is all the more surprising given that the 
geometric distribution should better accommodate the 
dititeouet in capture probability between males and 
females than should the Poisson distribution. 

The Lincoln-Petersen method requires capture on 
two or more, and Jolly-Seber on three or more 
occasions. Both provide good population estimates 
provided assumptions are met (Southwood 1994). In 
our study, it is unlikely that we meet either the 
assumptions of a closed population or equal capture 
probability required by the Lincoln-Peterson method. 
For the Lincoln Peterson method, the loss and gain of 
individuals after the initial marking period will result in 
overestimation of population size (Gall 1985). Jollv- 
Seber has the advantage of providing parameters for 
capture, survival, and recruitment, but requires more 
sampling occasions. 


however 
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TABLE 2. The relationship among population estimation methods. The dependent variable is in columns and the independent variable is in 
rows. The regression equation (SE) is on top and statistics for the relationship in the bottom of each cell. All relationships were significant (P 
< 0.001). A °° indicates that the intercept differs significantly from 0 (P < 0.05). 


Transect Craig’s Geometric Number of Lincoln- Jolly-Seber — 
captures Petersen 
Transect 1.17£0.36** 1.48+0.40 0.9740,29** 1.8940.34** 1.54+0.46** 
1.0540.13 1.12x1.45 0.80+£0.11 0.80x0.12 0.7720.16 
R/-0.80 R?—-0.79 R/-0.78 R/-0.81 R’=0.64 
F116764.2 F146759.7 F116757.3 F311746.41 F113722.7 
Craig's -0.4140.32 0.2040.04** 0.3140.33 0.954£0.30** 0.83+0.48 
0.77x0.10 1.0820.01 0.70+0.08 0.79+0.07 0.70+0.12 
R?-0.80 R?-0.99 R?-0.82 R?-0.91 R?20.73 
F116764.2 F11671093.2 F116773.5 F1117113.5 Fj137235.2 
Geometric -0.18x0.40** 0.1940.34 0.7840.32** 0.72+0.50 
0.93+0.01 0.65+£0.08 0.7423:0.07 0.6520.11 
R?-0.99 R?-0.82 R?-0.91 R?=0.73 
F11671093.2 F146772.0 Fi 115112.3 Fi 13=34.6 
Number of -0.41+0.41 0.3040.43 0.53x0.47 0.7440.31 0.52x0.34 
captures 0.980. 13 1.17+0.14 1.2640.15 1.04+0.10 1.0040.11 
R’=0.78 R’=0.82 R/-0.82 R?-0.92 R’=0.87 
JON ae os) ieee ion) F\ 16=72.0 F117118.7 F143785.8 
Lincoln- -1.42x0.62** -0.77x0.45 -0.58:0.48 -0.39x0.33 0.16x0.60 
Petersen 1.0140.15 1.16x0.11 1.24+0.12 0.88+0.08 0.86+0.14 
R/-0.81 UP R^-0.91 R?-0.92 R/-0.76 
F111746.4 Gmie SS) F1117112.3 Fi 115118.7 Fy „1=35.5 
Jolly- -0.3740.64 0.14+0.65 0.3440.70 -0.074£0.34 0.804£0.56 
Seber 0.83+0.17 1.05£0.18 1.13+0.19 0.87£0.09 0.89+40.15 
R/-0.64 R?-0.73 R’=0.73 R°=0.87 R/-0.76 
ATE di Fy 13=35.2 F\13=34.6 Fi 13=85.8 Fi. 111735.5 
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